Entanglement and boundary critical phenomena 
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We investigate boundary critical phenomena from a quantum information perspective. Bipartite 
entanglement in the ground state of one-dimensional quantum systems is quantified using the Renyi 
entropy S a , which includes the von Neumann entropy (a — > 1) and the single-copy entanglement 
(a — > oo) as special cases. We identify the contribution from the boundary entropy to the Renyi 
entropy, and show that there is an entanglement loss along boundary renormalization group (RG) 
flows. This property, which is intimately related to the Afneck-Ludwig p-theorem, can be regarded as 
a consequence of majorization relations between the spectra of the reduced density matrix along the 
boundary RG flows. We also point out that the bulk contribution to the single-copy entanglement 
is half of that to the von Neumann entropy, whereas the boundary contribution is the same. 
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Recently much work has been done to understand en- 
tanglement in quantum many-body systems. In partic- 
ular, the behavior of various entanglement measures at 
or near a quantumphase transition [jj has received a 

lot of attention |lllSIEIS0EI3- TheSe entan g le " 
ment measures include the von Neumann entropy and the 
single-copy entanglement, among others. The former is 
the most studied measure and quantifies entanglement in 
a bipartite system in the so-called asymptotic regime |lfjj , 
whereas the latter was recently suggested to quantify the 
entanglement present in a single copy 9]. For a system 
in a pure state (e.g. the ground state) that is parti- 
tioned into two subsystems A and B, the von Neumann 
entropy is Si = -TrAPAiog 2 pa where pA = Tr s |^)(i/i| 
is the reduced density matrix for A, and the single-copy 
entanglement is Soo = — log 2 Ai, where Ai is the largest 
eigenvalue of pa- 

Studies of the von Neumann entropy for quantum spin 
chains 0, S IS IE 13 have revealed that its depen- 
dence on the size £ of the block A is very different for 
noncritical and critical systems. For the former, the 
von Neumann entropy increases logarithmically with £ 
until it saturates when £ becomes of order the correla- 
tion length £, while for the latter (having £ = oo) it 
diverges logarithmically with I. Remarkably, the prefac- 
tor of the logarithmic term is universal and proportional 
to the central charge of the underlying conformal field 
theory (CFT) Id. 1 1 If. Furthermore, it has been shown 
numerically |l2| that the entanglement loss along the 
(bulk) renormalization group (RG) flows, which is con- 
sistent with the CFT predictions for the von Neumann 
entropy and with Zamolodchikov's c-theorem 

can be given a more "fine-grained" characterization in 
terms of the majorization concept [lij . A theoretical 
analysis of majorization in these systems also appeared 
recently (lif . 

Boundary critical phenomena [lfij ] in one-dimensional 



(ID) quantum systems (equivalently, 2D classical sys- 
tems) have attracted a lot of interest, especially in the 
context of boundary CFT. A closely related subject is the 
theory of boundary perturbations of certain conformally 
invariant theories, so-called integrable boundary quan- 
tum field theory |l7|. which is relevant to quantum spin 
chains with nontrivial boundary interactions, impurities 
in Luttinger liquids, Kondo physics, tunneling in frac- 
tional quantum Hall devices, and open string theory. In 
these problems, the Affleck-Ludwig (/-theorem P ia Y s 
a similar role as the c-theorem does for bulk critical phe- 
nomena. 

In this Letter we investigate boundary critical phenom- 
ena in quantum spin chains from a quantum information 
perspective. Our findings provide new insights into the 
information-theoretic explanation of the boundary en- 
tropy and the ^-theorem characterizing the intrinsic ir- 
reversibility due to information loss along boundary RG 
flows [19j. We choose the Renyi entropy (also known as 
the a-entropy) S a = (1 — a) -1 log 2 TtaPa as our entan- 
glement measure, partly motivated by the fact that both 
the von Neumann entropy and the single-copy entangle- 
ment are special cases of the Renyi entropy, correspond- 
ing to a — > 1 and a — > oo, respectively. Using CFT we 
derive expressions for the Renyi entropy which include 
a boundary entropy contribution, and show that a ma- 
jorization relation underlies an entanglement loss along 
boundary RG flows (i.e., when a system with a bound- 
ary interaction flows from an unstable to a stable fixed 
point). We support our analytical arguments with nu- 
merical density-matrix renormalization group (DMRG) 
calculations [2(il2ll|. 

Renyi entropy and conformal field theory. — Consider 
a ID lattice of interacting spins with lattice spacing a. 
Let L be the total length of the system, and let the two 
subsystems A and B be blocks of consecutive spins of 
length £ and L — £, respectively. Then CFT predicts that 
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for an infinite spin chain at criticality, 

S a = ^{l + a- 1 )log 2 -+c' a (f) 
6 a 

where c is the central charge and c' a is a non-universal 
constant [|J . For a semi-infinite spin chain with the block 
A starting at the origin where a boundary interaction is 
applied, we have instead 

r 9f 1 

S Q = -(l + a-')log 2 - + -c' a + S b . (2) 
12 a 2 

Here Sb = log 2 g is the boundary entrop y ll Si . w ith g = 
(B\0), where \B) is a boundary state [22|, |23|, |24| and |0) 
is the ground state. We emphasize that, compared with 
the corresponding expression for Si in Ref. |f|, there is 
an extra factor 1/2 in front of c' a and Sb in Eq. 

The corresponding Renyi entropies for finite L are 
found from Eqs. and (J2J by standard conformal map- 
pings 23]. As a result, the Renyi entropy for periodic 
boundary conditions (PBC) [open boundary conditions 
(OBC)] is given by replacing l/a — > (L/na) smn£/L in 
Eq. Q [Eq. 101] • For the OBC case we have assumed 
that the BC on the left and right ends are identical; oth- 
erwise one has to consider the complicating effects of a 
boundary condition changing operator. 

From these results one sees that the contribution of 
the bulk universal part to the single-copy entanglement 
5*00 is always half of that to the von Neumann entropy 
Si, thus extending the conclusion in Ref. for the XX 
chain (c = 1) to all conformally invariant critical systems. 
However, we stress that the contribution of the boundary 
entropy to S a does not depend on a. 

The details of the CFT derivation of the above results 
(which makes use of CFT expressions for Tr^p^ from 
Ref. |f|) will be given elsewhere. Here we will instead 
present a heuristic argument to justify Eqs. JTJ and J3J|. 
As is well known, the central charge c, which measures 
the (effective) number of gapless excitations, describes 
the way a specific system reacts to the introduction of 
a macroscopic length scale into the system [23j. There- 
fore, when the entire (infinite) system is partitioned into 
a block of length I and its environment, one may expect 
that the Renyi entropy S a depends only on c and l/a 
with some short distance cutoff a, e.g. the lattice spacing 
for quantum spin chains. Because both the Renyi entropy 
and the central charge are additive, one can see |25j that 
the Renyi entropy must be linear as a function of the cen- 
tral charge, i.e., S a = cf a {£/a) + h a , with f a a universal 
function and h a a non-universal (i.e., model-dependent) 
function. The specific form of the function f a may be de- 
termined by calculating the Renyi entropy for any exactly 
solvable model, for instance the massless Dirac fermion 
field as done in Ref. [2j for PBC, thus confirming Eq. @J. 
The calculation may be extended to the massless Dirac 
fermion with a conformally invariant boundary, leading 
to the conclusion that h a consists of the non-universal 



bulk part 1/2 c' a and the universal boundary part Sb- 
The bulk universal part of Eq. © and its counterpart 
for finite L result from the substitution c — > c/2, I — > 21 
and L — > 2L in Eq. Q and its corresponding counter- 
part, due to the fact that one may "unfold" the system 
by identifying left-movers at position x (x = ja, with j 
labeling lattice sites) with right-movers at —x so that the 
resulting system consists of only right movers subject to 
PBC and with the system size scales I and L doubled 
and the number of gapless degrees of freedom halved. 
The factor 1/2 in front of c' a in the second bulk term in 
Eq. (J2J also originates from this halving of the number 
of gapless degrees of freedom. 

S — 1/2 transverse Ising and XXZ chains in bound- 
ary magnetic fields. — We now address the implications of 
these CFT predictions for models of semi-infinite quan- 
tum spin chains. We first consider the S — 1/2 transverse 
Ising chain in a boundary magnetic field, described by 

oo 

-fusing = — y^(SJSJ + i + hSj) + hbS$. (3) 

3=0 

Here h is the transverse bulk magnetic field, and hb is 
the boundary magnetic field. We set h = 1/2 so that the 
model is bulk critical with central charge c = 1/2. The 
second model is the S — 1/2 XXZ chain, for which 

oo 

ffxxz - + s y s y +i + AS|S| + i) + h h S$. (4) 

j=o 

Here A denotes the anisotropy and hb is the (transverse) 
boundary magnetic field. The model is bulk critical with 
central charge c = 1 for — 1 < A < 1. 

For both models, the points hb = and hb = ±oo, 
corresponding to the free and fixed conformally invari- 
ant BC, are boundary critical fixed points; the former is 
unstable and the latter is stable. A boundary magnetic 
field hb > is generally a relevant perturbation of the 
free BC hb = 0, and generates a boundary RG flow of 
hb towards the fixed point hb — oo. For < hb < oo, 
the conformal invariance is lost, and the competition be- 
tween boundary ordering and bulk criticality introduces 
a crossover length ( cx h^~ x j2||, where d < 1 is the 
scaling dimension of the relevant boundary perturbation. 
For the Ising model, d — 1/2, and for the XXZ model, 
d = 2uR 2 , where R = ^{l/2ir) - (1/2tt 2 ) arccos A is the 
compactification radius. Furthermore, for both models, 
g for fixed BC is less than g for free BC, which implies 
that the Renyi entropy is less for fixed than for free BC. 
This is also consistent with the g-theorem 0], which 
states that g decreases along boundary RG flows. For 
the transverse Ising model, g = 1 (free) and g = l/\[2 
(fixed) 0; for the XXZ model, g = 7 r^ 1 / 4 (2i?)^ 1 / 2 (free) 
and g = ix 1 I^R}l' 1 (fixed) |2^|. We emphasize that even 
away from boundary critical points, Eq. @ is still valid 
for I > C, due to bulk criticality. 
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FIG. 1: Bipartite entanglement for the bulk critical transverse 
Ising model, quantified by the von Neumann entropy Si of the 
reduced density matrix. The block size £ is given as T(£) (see 
text) which allows for a direct check of the CFT predictions 
Q and (j2J. The fit yields the predicted boundary entropies. 



Hxxz reduces to the isotropic XXX (Heisenberg) 
model for A = 1 (R = \/\/2rv). This case is special, 
because the boundary perturbation is marginal (d = 1). 
In fact, g = 2 -1 / 4 for both free and fixed BC. The line 
from hi, = to hb — oo is a line of fixed points; there is 
no RG flow since g is the same everywhere along the line. 
To see this, we consider the finite XXX chain with iden- 
tical boundary interactions at the two boundaries. This 
corresponds to a free scalar boson field with a dynamical 
boundary interaction |2t1 ]. for which it was found that 
when both ends have the same BC, the partition func- 
tion is independent of the boundary interaction strength, 
reflecting the fact that the energy levels do not feel the 
presence of the boundary interaction. Combining this 
with the CFT prediction S = {nc/3f3)L + 2S b for 
the Gibbs entropy S = —f3 2 dF/d(3, we conclude that 
the boundary entropy does not depend on the boundary 
magnetic field hb- 

To check the CFT predictions, we analyze quasi-exact 
DMRG results. Fig. ^ shows the entanglement in the 
bulk critical transverse Ising model © for PBC and 
OBC with different boundary fields hi,. For PBC we use 
the DMRG variant described in Ref. j^. The block 
size £ is shown as T(£), with T PB cM = log 2 sin ^) 

and TobcW = \ log 2 (§§ sm T")' a PP ro P r i ate f° r a un ~ 
ear fit of the finite L counterparts of J3> and @. The 
fits for PBC, free OBC and fixed OBC yield the de- 
sired central charge (c = 0.500, 0.502, 0.498 » 1/2) and 
the predicted boundary entropies (S^ ee = —0.003 « 0, 
5 fixed = _ _ 4 g 7 ~ log 2 l/V2), which can be read off 
from the differences of the axis intercepts. The curves for 
nonzero hi, converge for large £ to the fixed BC curve (the 
smaller the value of hi,, the larger the crossover length £). 
This implies that the boundary entropy for any nonzero 



hi, is the same as that for fixed BC, consistent with the 
explicit construction of boundary states in Ref. ^3 ■ 

Majorization relations. — Let us first recall the defi- 
nition of majorization |l4|. Consider two probability 
distributions A = {Ai} and p, = {pi} whose elements 
are ordered such that Ai > A 2 > • • • > A r , and simi- 
larly for p. We say that A is majorized by /i, written 
A -< p,, if Yt,i=i - 2i=i Mi f° r k = 1, . . . , r - 1 and 
Y^i=i — Y^i=i ( = !)• Here the probability distribu- 
tions are formed by the eigenvalues of the reduced density 
matrix oa , and r is the Schmidt rank. An elementary re- 
sult [1J| states that A -< p, if and only if 0(A) > cj>{p) 
for all Schur-concave functions </> In fact, the Renyi 
entropy is Schur-concave for any index a. Therefore it is 
necessary (but not sufficient) for the Renyi entropy for 
all indices a to be monotonic with some system parame- 
ter (£ or hb) for the corresponding spectra of the reduced 
density matrix to be subject to a majorization relation. 

The CFT results (JTJ and © show that at (both bulk 
and boundary) criticality the Renyi entropy increases 
monotonically with the block size £ (for this to hold for a 
finite chain of length L, £ must be less than L/2). In par- 
ticular, the largest eigenvalue Ai = 2~ s °° of pa decreases 
with increasing £. This indicates that pa' -< Pa if the 
block A is a sub-block of A 1 31]. Indeed, the majoriza- 
tion relation follows from the fact that any two eigenvalue 
distributions, corresponding to two different block sizes, 
only cross once when the eigenvalues A.; are plotted ver- 
sus the eigenvalue index i. That is, an index i* exists, 
such that Ai decreases with increasing £ for i < i*, and 
Ai increases with increasing £ for i > i* . The unique- 
ness of this crossing is guaranteed by the fact that the 
eigenvalue distribution follows Ai oc q l3i for both bulk 
and boundary conformal theories, apart from degenera- 
cies rii. Here q — exp(27rzT) for conformal theories and 
q = cxp(— tv 5) for boundary conformal theories, and (3i 
are integers related to the scaling dimensions of the de- 
scendant operators. The conformal invariance requires 
that both t and 8 should be proportional to 1/ In £. The 
discreteness of and Pi ensures that they do not change 
when the block size £ is varied (for finite L, £ must be 
restricted to even or odd values if the model has a parity 
effect). 

Next, we consider the behavior of the Renyi entropy 
along boundary RG flows. By combining our CFT re- 
sults with the (/-theorem, it follows that the Renyi en- 
tropy decreases (more precisely, does not increase) along 
boundary RG flows. In particular, the largest eigenvalue 
Ai of pa increases along the boundary RG flow. Fur- 
thermore, one may argue that along the boundary RG 
flow, the eigenvalues of pa take the same form Ai a q^ 
as at the conformally invariant fixed point at the end of 
the flow, except that the dependence of q on £ is differ- 
ent (again, the discreteness of the degeneracies ni and 
Pi ensures that they remain the same along the flow). 
However, for any nonzero h b , the bulk criticality requires 
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FIG. 2: Eigenvalue distributions for the reduced density ma- 
trix of a block of size I — 128a in the bulk critical transverse 
Ising model Q for various values of the boundary magnetic 
field hb- The quasi-exact DMRG results confirm the majoriza- 
tion relation along the boundary field RG flow (see text) . The 
inset shows that the crossing of the spectra occurs at i* — 1. 



that the dependence S ~ 1/ \n£ is recovered for £ > C. It 
then follows that there is one and only one crossing for 
the eigenvalue distributions versus i along the boundary 
RG flows. More precisely, an index i* exists, such that 
Xi increases for i < i* , and decreases for i > i* along the 
flows. This in turn implies the majorization relation. 

A rough estimate for the crossing index i* can be ob- 
tained from the CFT predictions for the largest eigen- 
value Ai. This gives i* ~ £ c / 6 for PBC and i* ~ £ c ' 12 for 
OBC. Thus the block size £ must be sufficiently large to 
observe that the crossing occurs at i* > 1. For instance, 
to see that the second largest eigenvalue decreases with 
increasing £, £ should at least be 1800a for the semi- 
infinite XX chain, as estimated from the exact solution 
[32^| . DMRG calculations of the reduced density ma- 
trix spectra show majorization along boundary RG flows, 
where eigenvalues down to ~ 10~ 15 are considered. Fig. [3 
shows the crossing point i* of the spectra in the bulk crit- 
ical transverse Ising model with block size £ = 128a and 
various values of the boundary field h^. The crossing of 
the spectra occurs here at i* = 1. 

For the Heisenberg XXX model (A = 1) the Renyi en- 
tropy does not depend on the boundary magnetic field: 
dS a /dhh = 0. These constraints (of which there is an in- 
finite number, one for each a) imply that the spectra of 
Pa do not vary with h^. Therefore the presence of a line 
of fixed points here amounts to the statement that all of 
these critical points share the same eigenvalue distribu- 
tion, so there is no boundary RG flow (a boundary RG 
flow would require a preferential direction characterized 
by (strict) majorization). 

Conclusions. We have investigated the interrelations 
between the boundary entropy, the Affleck-Ludwig g- 



theorem and entanglement in quantum spin systems with 
boundary interactions. The intrinsic irreversibility along 
boundary RG flows, as embodied in Affleck-Ludwig g- 
theorem, is connected with the majorization relation 
solely characterized by the ground state itself. The re- 
sults will bring new insights into our understanding of 
2D statistical mechanical systems with boundary inter- 
actions and (perturbed) conformal theories defined on 
manifolds with boundaries. The equivalence between ID 
quantum systems and 2D classical systems indicates that 
majorization relations should still be preserved in 2D 
classical systems, even if it does not make sense to speak 
of quantum entanglement for classical cases. 

During the preparation of this manuscript, two 
preprints [3^] appeared which reached the same conclu- 
sion as us regarding the connection between the von Neu- 
mann entropy and the single-copy entanglement for bulk 
conformal theories. We thank Guifre Vidal for enlight- 
ening discussions and comments. HQZ and JOF thank 
the Australian Research Council for support. TB thanks 
the DFG for support. 
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